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Abstract 

We present the semileptonic decay of B c meson in a QCD potential 
model with the coulombic part of the Cornell potential — +br+c as 
perturbation. Computing the slope and curvature of Isgur Wise func- 
tion in this approach, we study the pseudoscalar and vector form fac- 
tors for the transition of B c meson to its S wave charmonium cc states. 
Numerical estimates of widths for the transitions of B c — > J /ip(ri c )lvi 
are presented. The results are found to be in good agreement with 
other theoretical values. 
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1 Introduction 

The investigation of weak decays of mesons composed of a heavy quark 
and antiquark gives a very important insight in the heavy quark dynam- 
ics. The exclusive semileptonic decay processes of heavy mesons generated a 
great excitement not only in extracting the most accurate values of Cabbibo- 
Kobayashi Maskawa(CKM) matrix elements but also in testing diverse the- 
oretical approaches to describe the internal structure of hadrons. The great 
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virtue of semileptonic decay processes is that the effects of the strong in- 
teraction can be separated from the effects of the weak interaction into a 
set of Lorentz- invariant form factors, i.e., the essential informations of the 
strongly interacting quark/gluon structure inside hadrons. Thus, the the- 
oretical problem associated with analyzing semileptonic decay processes is 
essentially that of calculating the weak form factors. 

The decay properties of thel? c meson are of special interest, since it is the 
only heavy meson consisting of two heavy quarks with different flavor. This 
difference of quark flavors forbids annihilation into gluons. As a result, the 
excited B c meson states lying below the Bp meson threshold undergo pionic 
or radiative transitions to the pseudoscalar ground state which is consider- 
ably more stable than corresponding charmonium or bottomonium states and 
decays only weakly. The CDF Collaboration reported the discovery of the B c 
ground state in pp collisions already more than ten years ago [1]. However, up 
till recently its mass was known with a very large error. Now it is measured 
with a good precision in the decay channel B c — > J/ipn. More experimental 
data on masses and decays of the Bc mesons are expected to come in near 
future from the Tevatron at Fermilab and the Large Hadron Collider (LHC) 
at CERN. The estimates of the Bc decay rates indicate that the c quark 
transitions give the dominant contribution while the b quark transitions and 
weak annihilation contribute less. However, from the experimental point of 
view the B c decays to charmonium are easier to identify. Indeed, CDF and 
DO observed the B c meson and measured its mass analyzing its semileptonic 
and nonleptonic decays B c — > J/iplu. 

There are many theoretical approaches to the calculation of exclusive B c 
semileptonic decay modes. Some of them are: QCD sum rules [ 2, 3, 4], the 
relativistic quark model [5, 6, 7] based on an effective Lagrangian describ- 
ing the coupling of hadrons to their constituent quarks, the quasipotential 
approach to the relativistic quark model [8, 9, 10], the instantaneous non- 
relativistic approach to the Bethe-Salpeter(BS) equation [11], the relativistic 
quark model based on the BS equation [12, 13], the QCD relativistic potential 
model [14], the relativistic quark- meson model [15], the nonrelativistic quark 
model [16], the covariant light-front quark model [17], and the constituent 
quark model [18, 19, 20, 21] using BSW(Bauer, Stech, and Wirbel) model 
[22] and ISGW(Isgur, Scora, Grinstein, and Wise) model [23]. The purpose 
of this paper is to extend a QCD potential model [24] with coulombic part as 
perturbation to calculate the hadronic form factors and decay widths for the 
exclusive semileptonic decay of B c meson. 
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Recently, we have reported the slope and curvature of I-W function for D 
and B mesons with the coulombic part of the potential as the perturbation 
in two particular renormalisation schemes MS and V scheme [24]. Instead 
of using a particular renormalisation schemes, here in this manuscript we use 
the strong coupling constant as a scale dependent parameter and compute 
the slope and curvature of I-W funtion for B c meson with a different set of 
mass input parameters than that of ref. [24]. We then use the I-W function 
to study the form factors and decay rates of B c meson into its 5* wave char- 
monium cc states within the framework of the potential model. 

The rest of the paper is organised as follows : section 2 contains the 
formalism with its subsections containing the model wavefunction, Masses, 
hadronic form factors and decay widths. In section 3 we place our results 
and conclusions. 



2 Formalism 

2.1 The wavefunction 

The wavefunction computed by Dalgarno method [25,26] with coulombic part 
— Tfr 1 + c °f the potential as perturbation and linear part br as parent has 
been reported in ref. [24]. For completeness we summarise the main equations 
in this section. However the alternate approach of choosing the linear part 
as perturbation has been reported earlier [27,28]. 

The total wave function corrected upto first order with normalisation is [24] 



where Ni is the normalisation constant for the total wave function ip cou i (r) 
with subscript 'coul' means coulombic potential as perturbation, 
where p 0n are given as [29[ [3U] : 
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and 



ba 4 x W 1 x 0.8808 x (bp)* 0.6535 x (b/i) s 
01 ~ (E-c) + 3a s (£ - c) (E^c) ( ' 

AnW 1 x 0.1183 
° 2 = 3^ (6) 

W« = / r 2 H' V (0) (r) rfr (7) 

JO 

*=-(£)'*• (8) 

here b and c are the model input parameters as is used in our previous works 
[24,27,28]. 'n' is the principal quantum no.(n=l for ground state), \x is the 
reduced mass of mesons and a s is the strong running coupling constant. 

2.2 Estimation of meson masses 

The energy shift of mass splitting due to spin interaction in the perturbation 
theory reads [31] 



leading to 



Taking the energy shift into account, meson mass is expressed as 

M p = m l + m 3 -- H<M0)| 2 (11) 

2>rriimj 

M v = m l + m, + p^\^(0)\ 2 (12) 
yrriimj 

since (Si.Sj) = —3/4 for pseudoscalar meson and (Si.Sj) = 1/4 for vector 
meson . 
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In the model approach, the wave function at the origin \ip (0) | develops sin- 
gularity. Therefore for computing \ip (0) | 2 we use the phenomenological law 
given by Igi and Ono[32] as 

|^(0)| 2 = fc/i L52 (13) 

here k is a proportionality constant. In eq.(l), the strong coupling constant 
connected to the potential is a function of the momentum as 

(11 - in (&) 

where rif is the number of flavour and /ii is the renormlistion scale related to 
the constituent quark mass and A is the QCD scale which is taken as 0.150 
GeV by fixing a s = 0.118 at the Z boson mass(91 GeV). We use the most 
common renormalisation scale as fii = with rif = 4 [32,33] and then 

evaluate a s for B c , r\ c and J/ip mesons. 

For numerical calculation, we use the model input parameters as, m c = 
1.55 GeV and nib = 4.79GeV and k = 0.05 to compute the masses of B c , r\ c 
and J/ip mesons and are shown in table 1. 



Table 1: Masses of cb and cc state 



mesons 


|V (0) | 2 in GeV 3 


values of a s 


M p>v in GeV 


B c 


0.061 


0.22 


6.274 [our work ] 
6.277 Expt[31! 
6.270 [H] 
6.302 [33] 


Vc 


0.034 


0.25 


3.070 [our work ] 
2.980 ExptpS] 
2.979 |8] 
3.088 [36] 


J/i> 


0.034 


0.25 


3.109 [our work] 
3.097 Expt[3l] 
3.096 [8] 
3.168 [36J 
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2.3 Pseudoscalar and Vector form factors 



In the case of the final cc states corresponds to the J = 0, r\ c states, as 
the matrix element of any axial current between the two pseudoscalar 
mesons vanishes, only vector current contributes. Unlike in the case of 
electromagnetic current of the charged pions, here the vector current = 
crf^b is not conserved as q^V 1 oc (m& — m c ) ^ 0. So the matrix element of 
the hadronic current, V 1 between the two J p = CT mesons is expressed in 
terms of two form factors f±{q 2 ) as 



(Vc(p')\V»\B c (p)) = f + (q 2 )(p + p), + f4q 2 )(p-p 



(15) 



Where q = p — p = ki + &2 is the four momentum transfer and f+{q 2 ) 
and f-{q 2 ) are the dimensionless weak transition form factors corresponds to 
B c — > Vc> which are functions of the invariant q 2 . Here q 2 varies within the 
range m 2 < q 2 < (m Bc - m Vc ) 2 = q 2 max . 

The transition between the pseudoscalar B c and the vector J/ip(p,e) 
mesons depends on four independant form factors as, 



{J/i/>(p,e)\crfb\B c (p)) 
(J/TP{p',e)\cj^ B b\B c (p)} 



M Bc + Mj,^ 



V(q 2 



(16) 



(M Bc + M 



-e • q 



(p + p 



e ■ qq h 



A,{q 2 



(M Bc - M JH ,)<f 



2M 



M Bc + M J/i3 
6 ■ qq^ 1 . / 2\ 

Mi ) 



A 2 (q 2 ) 
(17) 



In the present study we treat B c system similar to D (cd) system as in ref. [37] 
and extend HQET to express Eqns. [15] and [16] as 



M Bc M v< 



:(r lc (v')\V»\B c (v)) = (v + v'n(u;) 



1 



yjM Bc M JH , 



(J/^(v,e 3 )\V»\B c (v)) = ie^e v v' a v^(u) 



iM Bc Mj/^ 



(J/^v, e*)\A*\B c (v)) = [(1 + coy - (e ■ v)v'^(oj) 



(18) 

(19) 
(20) 
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Where £(cu) is the Isgur Wise function. For small, nonzero recoil, it is conven- 
tional to write the Isgur- Wise function as 



av.v') = ch 



i-p 2 (uj-i) + c(uj~iy + .. 



(21) 



where u is given by, 



w = v ■ v 



m l c + m lc ~ Q 2 



2m Bc m c 



[22) 



The quantity p 2 is the slope of I-W function at u — 1 and known as 
charge radius : 



2^1 



(23) 



The second order derivative is the curvature of the I-W function known as 
convexity parameter : 



(24) 



2 \du 2 / 

For the heavy-light flavor mesons the I-W function can also be written as 



£ (cu) = / 47rr |'0( r )| cosprdr 
Jo 



(25) 



where 



p 2 = 2/i 2 (w - 1) (26) 

Using eq.l,eq.21 and eq.25, we compute the slope and curvature of I-W func- 
tion as 

p\ = 0.67 and C Bc = 0.06. 
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Consequently, the form factors f±(q 2 ) correspond to the cc(rj c ) final state 
are related to the Isgur Wise function as [23, HI] 

/ ± tf)=«M^S (27) 



m 



Vc 



And those related to the J /if} as the final hadronic state are given by 

q2 rW*%^<«) 



V(q 2 ) = A 2 (q 2 ) = A (q 2 ) 



(M Bc + M J/V) ) 2 J 7 4M Bc M JV , 

(28) 

It is evident from Eqn. [251 that at q 2 — > 

ng 2 )=A 2 (g 2 )=A(g 2 )=A (g 2 ) (29) 

Thus knowing the masses and Isgur Wise function of the transition B c — > 
cc(r) c , J/if>)£ + isi, we will be able to compute respective form factors corre- 
spond to B c — > cc(r) c , J/if>)£ + U£ transitions. 

It should be noted that by virtue of transversality of the lepton current 
= ljfj,(l + 75)^ in the limit mi — > 0, the probabilities of semileptonic 
decays into e + v e and are independent of /_ . Thus, in calculation of 

these particular decay modes of B c meson this form factor can be consistently 
neglected [3]. 

The differential semileptonic decay rates can be expressed in terms of 
these form factors by 

(a) B c — > Pev decay (P = r] c ) 

AT C 2 A 3 IT/ I 2 

^_ (Bc _> Pev) = ^l^L\ S+{q ^. (30) 

(b) B c — > Vev decay (V = J/if>) The decays rate in transversally and longi- 
tudinally polarized vector mesons are defined by[42] 

dr L G 2 F A\V qb \ 2 q 2 2 2 

d?~ 96vr3 Ml l oiq)1, { ) 



dq 2 dg 2 dg 2 g 67r 3 M 



B 
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where helicity amplitudes are given by the following expressions 

2M Bc A 



M Bc + M v 



^%^ 2 > 



(33) 



H (q 2 



1 



2M yv V 

Thus the total decay rate is given by [42] 



a-\jt2 a 2 

(M Bc + M v )(M 2 Bc -Ml- q 2 )A x {q 2 ) - B 



M Bc + Mi 



v 



(34) 



96tt 3 M 



where G F is the Fermi constant, is CKM matrix element (q — c), 
A = |A 



(35) 



(Ml c + v - q*)* 



Integration over q 2 of these formulas gives the total rate of the corresponding 
semileptonic decay The computed results for the form factors and the total 
decay rate is shown in table 2, table 3 and table 4 respectively and compared 
with the other theoretical values. 



Table 2: Form factors for B c — > r] c £ + U£ 



work 


A(o) 


/-(0) 


/ +^Slmax) 


f (q 2 ) 


our work 


0.886 


-0.304 


1.065 


-0.365 


Vinodkumar[37] 


0.71 


-0.25 


0.78 


-0.28 


Faustov[42 ] 


0.47 




1.07 




Choi and Ji[43] 


0.546 




1.035 





3 Results and Discussion 

In this paper, we have computed the slope and curvature of Isgur Wise func- 
tion for B c meson, considering the coulombic part of the Cornell potential as 
perturbation. Presumably, the values of p 2 and C are found to be smaller 
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Table 3: Form factors for B c — > J/ip£ + i>i 



work 


V(0) 


^1 ^Slmax) 


Via 2 ) 


Our work 


0.891 


0.941 


1.0622 


Vinod kumar[37] 


0.70 


0.75 


0.77 


Faustov[42] 


0.49 


0.88 


1.34 



Table 4: Decay width for B c ->• cc(£ + u e ) in 10~ 1 5 GeV 
B c -> vd^i) B c -> J/t/j(£+u t ) 





16.4 


39.8 


[6] 


10.7 


28.2 




5.9 


17.7 


m] 


14.2 


34.4 


05] 


11.1 


30.2 


[35] 


8.31 


20.3 


US] 


11 ± 1 


28 ±5 


[47] 


2.1 


21.6 




8.6 


17.2 


1491 


10 


42 
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than the other theoretical values for D and B mesons. The reason is sup- 
posed to be due to the truncating of higher terms in the Airy function(which 
is an infinite series). To avoid the singularity in the wavefunction at the 
origin, we also use a model independent phenomenological parametrisation 
with the same mass input parameters to calculate the masses of cb and cc 
states in the ground state and is shown in table. 1. These masses are in good 
agreement with the experimental result. 

The computed I-W function and masses are used to study the hadronic form 
factors using HQET in analogy to D meson system and is compared with 
other available results which are shown in table. 2 and table. 3. The decay 
widths for B c — > cc{( + v^) and B c — > J/ip(£ + ue) are shown in table. 4 and 
found to be in agreement with other results. 
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